Broadly Applicable Control Approaches Improve Accuracy of ChIP-Seq Data.

Chromatin ImmunoPrecipitation (ChIP) is a widely used method for the analysis of protein-DNA interactions in vivo; however, ChIP has pitfalls, particularly false-positive signal enrichment that permeates the data. We have developed a new approach to control for non-specific enrichment in ChIP that involves the expression of a non-genome-binding protein targeted in the IP alongside the experimental target protein due to the sharing of epitope tags. ChIP of the protein provides a "sensor" for non-specific enrichment that can be used for the normalization of the experimental data, thereby correcting for non-specific signals and improving data quality as validated against known binding sites for several proteins that we tested, including Fkh1, Orc1, Mcm4, and Sir2. We also tested a DNA-binding mutant approach and showed that, when feasible, ChIP of a site-specific DNA-binding mutant of the target protein is likely an ideal control. These methods vastly improve our ChIP-seq results in S. cerevisiae and should be applicable in other systems.

DNA binding specificity of all four Saccharomyces cerevisiae forkhead transcription factors.

Quantifying the nucleotide preferences of DNA binding proteins is essential to understanding how transcription factors (TFs) interact with their targets in the genome. High-throughput in vitro binding assays have been used to identify the inherent DNA binding preferences of TFs in a controlled environment isolated from confounding factors such as genome accessibility, DNA methylation, and TF binding cooperativity. Unfortunately, many of the most common approaches for measuring binding preferences are not sensitive enough for the study of moderate-to-low affinity binding sites, and are unable to detect small-scale differences between closely related homologs. The Forkhead box (FOX) family of TFs is known to play a crucial role in regulating a variety of key processes from proliferation and development to tumor suppression and aging. By using the high-sequencing depth SELEX-seq approach to study all four FOX homologs in Saccharomyces cerevisiae, we have been able to precisely quantify the contribution and importance of nucleotide positions all along an extended binding site. Essential to this process was the alignment of our SELEX-seq reads to a set of candidate core sequences determined using a recently developed tool for the alignment of enriched k-mers and a newly developed approach for the reprioritization of candidate cores.

Dbf4 Zn-Finger Motif Is Specifically Required for Stimulation of Ctf19-Activated Origins in Saccharomyces cerevisiae.

Eukaryotic genomes are replicated in spatiotemporal patterns that are stereotypical for individual genomes and developmental profiles. In the model system Saccharomyces cerevisiae, two primary mechanisms determine the preferential activation of replication origins during early S phase, thereby largely defining the consequent replication profiles of these cells. Both mechanisms are thought to act through specific recruitment of a rate-limiting initiation factor, Dbf4-dependent kinase (DDK), to a subset of licensed replication origins. Fkh1/2 is responsible for stimulation of most early-firing origins, except for centromere (CEN)-proximal origins that recruit DDK via the kinetochore protein Ctf19, which is required for their early firing. The C-terminus of Dbf4 has been implicated in its recruitment to origins via both the Fkh1/2 and Ctf19 mechanisms. Here, we show that the Zn-finger motif within the C-terminus is specifically required for Dbf4 recruitment to CENs to stimulate CEN-proximal/Ctf19-dependent origins, whereas stimulation of origins via the Fkh1/2 pathway remains largely intact. These findings re-open the question of exactly how Fkh1/2 and DDK act together to stimulate replication origin initiation.

Rpd3 regulates single-copy origins independently of the rDNA array by opposing Fkh1-mediated origin stimulation.

Eukaryotic chromosomes are organized into structural and functional domains with characteristic replication timings, which are thought to contribute to epigenetic programming and genome stability. Differential replication timing results from epigenetic mechanisms that positively and negatively regulate the competition for limiting replication initiation factors. Histone deacetylase Sir2 negatively regulates initiation of the multicopy (∼150) rDNA origins, while Rpd3 histone deacetylase negatively regulates firing of single-copy origins. However, Rpd3's effect on single-copy origins might derive indirectly from a positive function for Rpd3 in rDNA origin firing shifting the competitive balance. Our quantitative experiments support the idea that origins compete for limiting factors; however, our results show that Rpd3's effect on single-copy origin is independent of rDNA copy-number and of Sir2's effects on rDNA origin firing. Whereas RPD3 deletion and SIR2 deletion alter the early S phase dynamics of single-copy and rDNA origin firings in opposite fashion, unexpectedly only RPD3 deletion suppresses overall rDNA origin efficiency across S phase. Increased origin activation in rpd3Δ requires Fkh1/2, suggesting that Rpd3 opposes Fkh1/2-origin stimulation, which involves recruitment of Dbf4-dependent kinase (DDK). Indeed, Fkh1 binding increases at Rpd3-regulated origins in rpd3Δ cells in G1, supporting a mechanism whereby Rpd3 influences initiation timing of single-copy origins directly through modulation of Fkh1-origin binding. Genetic suppression of a DBF4 hypomorphic mutation by RPD3 deletion further supports the conclusion that Rpd3 impedes DDK recruitment by Fkh1, revealing a mechanism of Rpd3 in origin regulation.

Dynamic relocalization of replication origins by Fkh1 requires execution of DDK function and Cdc45 loading at origins.

Chromosomal DNA elements are organized into spatial domains within the eukaryotic nucleus. Sites undergoing DNA replication, high-level transcription, and repair of double-strand breaks coalesce into foci, although the significance and mechanisms giving rise to these dynamic structures are poorly understood. In S. cerevisiae, replication origins occupy characteristic subnuclear localizations that anticipate their initiation timing during S phase. Here, we link localization of replication origins in G1 phase with Fkh1 activity, which is required for their early replication timing. Using a Fkh1-dependent origin relocalization assay, we determine that execution of Dbf4-dependent kinase function, including Cdc45 loading, results in dynamic relocalization of a replication origin from the nuclear periphery to the interior in G1 phase. Origin mobility increases substantially with Fkh1-driven relocalization. These findings provide novel molecular insight into the mechanisms that govern dynamics and spatial organization of DNA replication origins and possibly other functional DNA elements.

Quantitative Bromodeoxyuridine Immunoprecipitation Analyzed by High-Throughput Sequencing (qBrdU-Seq or QBU).

Incorporation into DNA of nucleoside analogs like 5-bromo-2'-deoxyuridine (BrdU) is a powerful tool for in vivo studies of DNA synthesis during replication and repair. Immunoprecipitation of BrdU-labeled DNA analyzed by DNA sequencing (BrdU-IP-seq) allows for genome-wide, sequence-specific tracking of replication origin and replication fork dynamics under different conditions, such as DNA damage and replication stress, and in mutant strains. We have recently developed a quantitative method for BrdU-IP-seq (qBrdU-seq) involving DNA barcoding to enable quantitative analysis of multiple experimental samples subjected to BrdU-IP-seq. After initial barcoding of multiple, individually BrdU-labeled genomic DNA samples, a pooling strategy is used for all subsequent steps including immunoprecipitation, amplification, and sequencing, which eliminates sample-to-sample variability in these steps. Parallel processing of an aliquot of the pooled input sample provides a direct control for the normalization of the data and yields results that allow quantitative comparisons of the experimental samples. Though developed for the analysis of S. cerevisiae, this method should be directly adaptable to other model systems.

Identification of Fkh1 and Fkh2 binding site variants associated with dynamically bound DNA elements including replication origins.

Forkhead Box (Fox) DNA binding proteins control multiple genome activities, including transcription, replication, and repair. These activities are organized spatially and temporally in the nucleus, and Fox proteins Fkh1 and Fkh2 have emerged as regulators of long-range chromosomal interactions involved with these activities, such as the clustering of replication origins programmed for early initiation. Fkh1 and Fkh2 bind a subset of replication origins and are thought to dimerize to mediate long-range chromosomal contacts between these origins. The binding of Fkh1 and/or Fkh2 (Fkh1/2) to replication origins and the recombination enhancer (RE), which is involved in DNA repair required for mating-type switching, is cell cycle-regulated and thus appears to be more dynamic than Fkh1/2 binding at regulated target genes. Here we report the identification of Fkh1/2 binding sequence variants at replication origins and the RE compared with Fkh1/2 binding sequences found at target genes of the CLB2 group. These different binding sequences have previously been characterized as weak and strong, respectively, suggesting that the presence of weak sites contributes to more dynamic interactions at replication origins and RE, possibly facilitated by Fkh1/2 dimerization and cooperative interactions with accessory proteins. We discuss the wealth of regulatory potential imbued in these features of the DNA and its binding proteins.

Reconstructing cell cycle pseudo time-series via single-cell transcriptome data.

Single-cell mRNA sequencing, which permits whole transcriptional profiling of individual cells, has been widely applied to study growth and development of tissues and tumors. Resolving cell cycle for such groups of cells is significant, but may not be adequately achieved by commonly used approaches. Here we develop a traveling salesman problem and hidden Markov model-based computational method named reCAT, to recover cell cycle along time for unsynchronized single-cell transcriptome data. We independently test reCAT for accuracy and reliability using several data sets. We find that cell cycle genes cluster into two major waves of expression, which correspond to the two well-known checkpoints, G1 and G2. Moreover, we leverage reCAT to exhibit methylation variation along the recovered cell cycle. Thus, reCAT shows the potential to elucidate diverse profiles of cell cycle, as well as other cyclic or circadian processes (e.g., in liver), on single-cell resolution.In single-cell RNA sequencing data of heterogeneous cell populations, cell cycle stage of individual cells would often be informative. Here, the authors introduce a computational model to reconstruct a pseudo-time series from single cell transcriptome data, identify the cell cycle stages, identify candidate cell cycle-regulated genes and recover the methylome changes during the cell cycle.

Conserved forkhead dimerization motif controls DNA replication timing and spatial organization of chromosomes in S. cerevisiae.

Forkhead Box (Fox) proteins share the Forkhead domain, a winged-helix DNA binding module, which is conserved among eukaryotes from yeast to humans. These sequence-specific DNA binding proteins have been primarily characterized as transcription factors regulating diverse cellular processes from cell cycle control to developmental fate, deregulation of which contributes to developmental defects, cancer, and aging. We recently identified Saccharomyces cerevisiae Forkhead 1 (Fkh1) and Forkhead 2 (Fkh2) as required for the clustering of a subset of replication origins in G1 phase and for the early initiation of these origins in the ensuing S phase, suggesting a mechanistic role linking the spatial organization of the origins and their activity. Here, we show that Fkh1 and Fkh2 share a unique structural feature of human FoxP proteins that enables FoxP2 and FoxP3 to form domain-swapped dimers capable of bridging two DNA molecules in vitro. Accordingly, Fkh1 self-associates in vitro and in vivo in a manner dependent on the conserved domain-swapping region, strongly suggestive of homodimer formation. Fkh1- and Fkh2-domain-swap-minus (dsm) mutations are functional as transcription factors yet are defective in replication origin timing control. Fkh1-dsm binds replication origins in vivo but fails to cluster them, supporting the conclusion that Fkh1 and Fkh2 dimers perform a structural role in the spatial organization of chromosomal elements with functional importance.

Quantitative BrdU immunoprecipitation method demonstrates that Fkh1 and Fkh2 are rate-limiting activators of replication origins that reprogram replication timing in G1 phase.

The Saccharomyces cerevisiae Forkhead Box (FOX) proteins, Fkh1 and Fkh2, regulate diverse cellular processes including transcription, long-range DNA interactions during homologous recombination, and replication origin timing and long-range origin clustering. We hypothesized that, as stimulators of early origin activation, Fkh1 and Fkh2 abundance limits the rate of origin activation genome-wide. Existing methods, however, are not well-suited to quantitative, genome-wide measurements of origin firing between strains and conditions. To overcome this limitation, we developed qBrdU-seq, a quantitative method for BrdU incorporation analysis of replication dynamics, and applied it to show that overexpression of Fkh1 and Fkh2 advances the initiation timing of many origins throughout the genome resulting in a higher total level of origin initiations in early S phase. The higher initiation rate is accompanied by slower replication fork progression, thereby maintaining a normal length of S phase without causing detectable Rad53 checkpoint kinase activation. The advancement of origin firing time, including that of origins in heterochromatic domains, was established in late G1 phase, indicating that origin timing can be reset subsequently to origin licensing. These results provide novel insights into the mechanisms of origin timing regulation by identifying Fkh1 and Fkh2 as rate-limiting factors for origin firing that determine the ability of replication origins to accrue limiting factors and have the potential to reprogram replication timing late in G1 phase.

ChIP-Seq to Analyze the Binding of Replication Proteins to Chromatin.

Chromatin immunoprecipitation (ChIP) is a widely used method to study interactions between proteins and discrete chromosomal loci in vivo. ChIP was originally developed for in vivo analysis of protein associations with candidate DNA sequences known or suspected to bind the protein of interest. The advent of DNA microarrays enabled the unbiased, genome-scale identification of all DNA sequences enriched by ChIP, providing a genomic map of a protein's chromatin binding. This method, termed ChIP-chip, is broadly applicable and has been particularly valuable in DNA replication studies to map potential replication origins in Saccharomyces cerevisiae and other organisms based on the specific association of certain replication proteins with these chromosomal elements, which are distributed throughout the genome. More recently, high-throughput sequencing (HTS) technologies have replaced microarrays as the preferred method for genomic analysis of ChIP experiments, and this combination is termed ChIP-Seq. We present a detailed ChIP-Seq protocol for S. cerevisiae that can be adapted for different HTS platforms and for different organisms. We also outline general schemes for data analysis; however, HTS data analyses usually must be tailored specifically for individual studies, depending on the experimental design, data characteristics, and the genome being analyzed.

Two-dimensional agarose gel electrophoresis for analysis of DNA replication.

The initiation, elongation, and termination of DNA replication are each associated with distinct, nonlinear DNA structures that can be resolved and identified by two-dimensional (2D) agarose gel electrophoresis. This method involves: isolation of genomic DNA while preserving fragile replication structures, digestion of the DNA with a restriction enzyme, separation of DNA by size and shape through two distinct stages of agarose gel electrophoresis, and Southern blotting to probe for the specific sequence(s) of interest. The method has been most commonly used to determine the activity level of putative replication origin-containing sequences, and has also been used to analyze replication timing, fork progression, fork pausing, fork stalling and collapse, termination, and recombinational repair.

Rif1 regulates initiation timing of late replication origins throughout the S. cerevisiae genome.

Chromosomal DNA replication involves the coordinated activity of hundreds to thousands of replication origins. Individual replication origins are subject to epigenetic regulation of their activity during S-phase, resulting in differential efficiencies and timings of replication initiation during S-phase. This regulation is thought to involve chromatin structure and organization into timing domains with differential ability to recruit limiting replication factors. Rif1 has recently been identified as a genome-wide regulator of replication timing in fission yeast and in mammalian cells. However, previous studies in budding yeast have suggested that Rif1's role in controlling replication timing may be limited to subtelomeric domains and derives from its established role in telomere length regulation. We have analyzed replication timing by analyzing BrdU incorporation genome-wide, and report that Rif1 regulates the timing of late/dormant replication origins throughout the S. cerevisiae genome. Analysis of pfa4Δ cells, which are defective in palmitoylation and membrane association of Rif1, suggests that replication timing regulation by Rif1 is independent of its role in localizing telomeres to the nuclear periphery. Intra-S checkpoint signaling is intact in rif1Δ cells, and checkpoint-defective mec1Δ cells do not comparably deregulate replication timing, together indicating that Rif1 regulates replication timing through a mechanism independent of this checkpoint. Our results indicate that the Rif1 mechanism regulates origin timing irrespective of proximity to a chromosome end, and suggest instead that telomere sequences merely provide abundant binding sites for proteins that recruit Rif1. Still, the abundance of Rif1 binding in telomeric domains may facilitate Rif1-mediated repression of non-telomeric origins that are more distal from centromeres.

Fkh1 and Fkh2 bind multiple chromosomal elements in the S. cerevisiae genome with distinct specificities and cell cycle dynamics.

Forkhead box (FOX) transcription factors regulate a wide variety of cellular functions in higher eukaryotes, including cell cycle control and developmental regulation. In Saccharomyces cerevisiae, Forkhead proteins Fkh1 and Fkh2 perform analogous functions, regulating genes involved in cell cycle control, while also regulating mating-type silencing and switching involved in gamete development. Recently, we revealed a novel role for Fkh1 and Fkh2 in the regulation of replication origin initiation timing, which, like donor preference in mating-type switching, appears to involve long-range chromosomal interactions, suggesting roles for Fkh1 and Fkh2 in chromatin architecture and organization. To elucidate how Fkh1 and Fkh2 regulate their target DNA elements and potentially regulate the spatial organization of the genome, we undertook a genome-wide analysis of Fkh1 and Fkh2 chromatin binding by ChIP-chip using tiling DNA microarrays. Our results confirm and extend previous findings showing that Fkh1 and Fkh2 control the expression of cell cycle-regulated genes. In addition, the data reveal hundreds of novel loci that bind Fkh1 only and exhibit a distinct chromatin structure from loci that bind both Fkh1 and Fkh2. The findings also show that Fkh1 plays the predominant role in the regulation of a subset of replication origins that initiate replication early, and that Fkh1/2 binding to these loci is cell cycle-regulated. Finally, we demonstrate that Fkh1 and Fkh2 bind proximally to a variety of genetic elements, including centromeres and Pol III-transcribed snoRNAs and tRNAs, greatly expanding their potential repertoire of functional targets, consistent with their recently suggested role in mediating the spatial organization of the genome.

The level of origin firing inversely affects the rate of replication fork progression.

DNA damage slows DNA synthesis at replication forks; however, the mechanisms remain unclear. Cdc7 kinase is required for replication origin activation, is a target of the intra-S checkpoint, and is implicated in the response to replication fork stress. Remarkably, we found that replication forks proceed more rapidly in cells lacking Cdc7 function than in wild-type cells. We traced this effect to reduced origin firing, which results in fewer replication forks and a consequent decrease in Rad53 checkpoint signaling. Depletion of Orc1, which acts in origin firing differently than Cdc7, had similar effects as Cdc7 depletion, consistent with decreased origin firing being the source of these defects. In contrast, mec1-100 cells, which initiate excess origins and also are deficient in checkpoint activation, showed slower fork progression, suggesting the number of active forks influences their rate, perhaps as a result of competition for limiting factors.

Location, location, location: it's all in the timing for replication origins.

The differential replication timing of eukaryotic replication origins has long been linked with epigenetic regulation of gene expression and more recently with genome stability and mutation rates; however, the mechanism has remained obscure. Recent studies have shed new light by identifying novel factors that determine origin timing in yeasts and mammalian cells and implicate the spatial organization of origins within nuclear territories in the mechanism. These new insights, along with recent findings that several initiation factors are limiting relative to licensed origins, support and shape an emerging model for replication timing control. The mechanisms that control the spatial organization of replication origins have potential impacts for genome regulation beyond replication.

Forkhead transcription factors establish origin timing and long-range clustering in S. cerevisiae.

The replication of eukaryotic chromosomes is organized temporally and spatially within the nucleus through epigenetic regulation of replication origin function. The characteristic initiation timing of specific origins is thought to reflect their chromatin environment or sub-nuclear positioning, however the mechanism remains obscure. Here we show that the yeast Forkhead transcription factors, Fkh1 and Fkh2, are global determinants of replication origin timing. Forkhead regulation of origin timing is independent of local levels or changes of transcription. Instead, we show that Fkh1 and Fkh2 are required for the clustering of early origins and their association with the key initiation factor Cdc45 in G1 phase, suggesting that Fkh1 and Fkh2 selectively recruit origins to emergent replication factories. Fkh1 and Fkh2 bind Fkh-activated origins, and interact physically with ORC, providing a plausible mechanism to cluster origins. These findings add a new dimension to our understanding of the epigenetic basis for differential origin regulation and its connection to chromosomal domain organization.

Strategies for analyzing highly enriched IP-chip datasets.

BACKGROUND: Chromatin immunoprecipitation on tiling arrays (ChIP-chip) has been employed to examine features such as protein binding and histone modifications on a genome-wide scale in a variety of cell types. Array data from the latter studies typically have a high proportion of enriched probes whose signals vary considerably (due to heterogeneity in the cell population), and this makes their normalization and downstream analysis difficult. RESULTS: Here we present strategies for analyzing such experiments, focusing our discussion on the analysis of Bromodeoxyruridine (BrdU) immunoprecipitation on tiling array (BrdU-IP-chip) datasets. BrdU-IP-chip experiments map large, recently replicated genomic regions and have similar characteristics to histone modification/location data. To prepare such data for downstream analysis we employ a dynamic programming algorithm that identifies a set of putative unenriched probes, which we use for both within-array and between-array normalization. We also introduce a second dynamic programming algorithm that incorporates a priori knowledge to identify and quantify positive signals in these datasets. CONCLUSION: Highly enriched IP-chip datasets are often difficult to analyze with traditional array normalization and analysis strategies. Here we present and test a set of analytical tools for their normalization and quantification that allows for accurate identification and analysis of enriched regions.